function [output] = steady_state_simulation_export_moments(parameters)
place = {'US', 'OECD'};
set_parameters
initialize_vectors

year = 1;

% MAIN SIMULATION LOOP
while year<=max_year;
    
    year_index = ['y' num2str(year)];
    if year > 1
        lag_year_index=['y' num2str(year-1)];
    end
    
    % SIMULATE INNOVATION PROCESS AND REASSIGN VARIETIES TO PRODUCERS
    innovation_and_assignment
        
    % ITERATE TO DETERMINE EQUILIBRIUM WAGES AND TRADE FLOWS
    [wage.OECD,omega,exports.US.(year_index),exports.OECD.(year_index),produced.US.(year_index),produced.OECD.(year_index),...
        nontraded.US.(year_index),nontraded.OECD.(year_index),prices.US.(year_index),prices.OECD.(year_index),output.US.(year_index),output.OECD.(year_index),...
        markups.US.(year_index),markups.OECD.(year_index),mu_bar.US(end+1),mu_bar.OECD(end+1),...
        nominal_GDP.US.(year_index),nominal_GDP.OECD.(year_index),export_share.US(end+1),export_share.OECD(end+1),revenue.US.(year_index),revenue.OECD.(year_index),x,m,...
        qualities_vec_consumed.US,qualities_vec_consumed.OECD,labor_share.US(end+1),labor_share.OECD(end+1),wage_level.US,wage_level.OECD,...
        P_bar.US.(year_index),P_bar.OECD.(year_index),total_income.US.(year_index),total_income.OECD.(year_index),...
        real_consumption.US.(year_index),real_consumption.OECD.(year_index),real_product_wage.US,real_product_wage.OECD] = ...
        find_equilibrium(wage.OECD,qualities_vec.US,qualities_vec.OECD,domestic_competitor_qual_vec.US,domestic_competitor_qual_vec.OECD,...
        import_tariff.US,import_tariff.OECD,l.US,l.OECD,n_varieties,export_share.US(end),export_share.OECD(end),mu_bar.US(end),mu_bar.OECD(end));
  
    if year == max_year-1
        calculate_trade_elasticity_real
    end
    
    % RUN FILE TO CALCULATE MOMENTS
    track_moments_within_loop    
       
    % UPDATE NEXT YEAR'S VARIETY OWNERSHIP BEFORE THE INNOVATION FOR THAT
    % YEAR BEGINS
    nxt_yr_index = ['y' num2str(year+1)];
    for cc=1:length(place)
    varown_vec.(place{cc}).(nxt_yr_index) = varown_vec.(place{cc}).(year_index);
    end
        
    l.US = population_growth*l.US;
    l.OECD = population_growth*l.OECD;
    
    year = year+1;
    
end

calculate_moments_after_loop

% CHOOSE A FEW KEY MOMENTS TO OUTPUT FROM THE FUNCTION

output = [mean_export_sales_growth_contracting_exporters_5yr_avg.US mean_export_sales_growth_contracting_exporters_5yr_avg.OECD; ...
mean_domestic_sales_growth_contracting_exporters_5yr_avg.US mean_domestic_sales_growth_contracting_exporters_5yr_avg.OECD; ...
mean_export_sales_growth_expanding_exporters_5yr_avg.US mean_export_sales_growth_expanding_exporters_5yr_avg.OECD; ...
mean_domestic_sales_growth_expanding_exporters_5yr_avg.US mean_domestic_sales_growth_expanding_exporters_5yr_avg.OECD];

end